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Abstract — In Cognitive Radio Networks (CRNs), unlicensed 
users are allowed to access the licensed spectrum when it is 
not currently being used by primary users (PUs). In this paper, 
we study the throughput maximization problem for a multi- 
channel CRN where each SU can only sense a limited number 
of channels. We show that this problem is strongly NP-hard, and 
propose an approximation algorithm with a factor at least |/j 
where fj, 6 [1,2] is a system parameter reflecting the sensing 
capability of SUs across channels and their sensing budgets. 
This performance guarantee is achieved by exploiting a nice 
structural property of the objective function and constructing 
a particular matching. Our numerical results demonstrate the 
advantage of our algorithm compared with both a random and 
a greedy sensing assignment algorithm. 

I. Introduction 

In the past decade, cognitive radio networks (CRNs) 
have emerged as a promising solution for achieving better 
utilization of the frequency spectrum to satisfy the increasing 
demand of wireless communication resources. In CRNs, 
secondary users (SUs) are offered the opportunity of access- 
ing the licensed channel when their activities do not cause 
disruptions for primary user (PU) transmissions. To this end, 
the Federal Communications Commission (FCC) [4] has 
opened the broadcast TV frequency bands for unlicensed 
users such as WLAN and WiFi. Most recently, congressional 
negotiators have reached the compromise to allow the auction 
of TV broadcast spectrum to wireless Internet providers [13]. 
IEEE has announced the IEEE 802.22 wireless network stan- 
dard [12] that specifies how to utilize the unused resources 
between channels in the TV frequency spectrum. 

To guarantee a high system throughput in a CRN, the main 
challenge is for the SUs to accurately detect the channel state 
of PUs while exploiting transmission opportunities over the 
white space. Sensing inaccuracies may lead to either a false 
alarm, where a channel is detected to be occupied when it is 
actually idle, or a misdetection, where a channel is detected 
to be idle when it is actually occupied. While the former hurts 
SU throughput, the latter hurts both PU and SU throughput. 
To improve sensing accuracy, cooperative spectrum sensing 
schemes [6], [9], [10] have been recently developed, where a 
joint decision is derived from individual observations made 
by multiple SUs, which effectively alleviates the impact of 
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incorrect individual decisions on throughput by exploiting 
the spatial diversity of the SUs. 

While cooperative sensing improves sensing accuracy, it 
also incurs sensing and reporting overhead at the SU side, 
especially when an SU senses multiple channels in a multi- 
channel CRN. In particular, requiring each SU to sense all 
the channels in a CRN may lead to long sensing durations, 
especially when the number of channels is large, which in 
turn reduces the average throughput of SUs. It is therefore 
reasonable to put a limit on the maximum sensing duration 
that an SU can afford, which translates to a budget on the 
number of channels that an SU can sense. Due to the hard- 
ware constraints, this budget could be different for different 
SUs. In this paper, we study the throughput optimization 
problem for a multi-channel CRN subject to this sensing 
constraint. 

Various cooperative sensing protocols have been proposed 
for maximizing system-wide performance metrics such as 
sensing accuracy [9] and system throughput [8], [16]. How- 
ever, these works either focus on a single-channel setting [9], 
[8] or allow each SU to sense all the channels [7], [3], [16]. 
In particular, an optimal Bayesian decision rule that maps a 
vector of local binary decisions made at SUs to a global 
decision on PU activity has been found for maximizing 
system throughput in a single channel setting [8], which 
achieves significantly better performance than linear rules 
such as AND, OR, and majority rules. However, a direct 
extension of the result in [8] to the multi-channel setting 
would require each SU to sense all the channels and incur 
high sensing duration. On the other hand, most works on 
multi-channel cooperative sensing put no explicit constraint 
on sensing duration of SUs. Furthermore, these works ei- 
ther use a simple linear decision rule [16] or require the 
transmission of the entire local sensing samples or sensing 
statistics at each SU. In our work, we choose to use a binary 
decision rule to avoid the high overhead involved in reporting 
complete local sensing results. However, instead of using a 
suboptimal linear rule as in [16], we use the optimal decision 
rule proposed in [8] for each channel. 

In this paper, we study the problem of maximizing the 
system throughput in a multi-channel CRN, by deciding 
for each channel, a subset of SUs to sense the channel, 
subject to the sensing budget constraint at each SU. Our 
main contributions can be summarized as follows: 

> We show that the throughput maximization problem is 
NP-hard in the strong sense and hence does not have a 
pseudo-polynomial time algorithm unless P = NP. 

« We prove that the system throughput function satisfies 



a structural property, and based on this we propose a 
matching-based algorithm, which achieves an approxi- 
mation factor at least \\i where /i £ [1,2] is a system 
parameter depending on the sensing capability of SUs 
across channels and their sensing budgets. 
This paper is organized as follows. The system model 
and the problem formulation are introduced in Section [II] 
In Section [HI] we prove that the optimization problem is 
NP-hard in the strong sense. We then prove the structural 
property of the system throughput function, and propose a 
matching based algorithm in Section [IV] In Section [V] nu- 
merical results illustrate the performance of our algorithms. 
The paper is concluded in Section [Vl] 

II. System Model 

In this section, we present the system model in two parts: 
communication model and cooperative sensing model. Based 
on the models, we formulate our overall objective, which is 
to decide the channel sensing assignment to maximize the 
overall system throughput. 

A. Communication Model 

We consider a time-slotted cognitive radio network com- 
posed of M orthogonal channels (each corresponding to a 
and N SUs. An SU may sense multiple PUs depending 
on its location. When the channel is idle, SUs that do 
not interfere with each other can transmit over it. Since 
scheduling and channel assignment for SU transmission are 
not the focus of this paper, we employ a simple policy: an 
SU is randomly selected for transmission over each available 
channel. Our model can readily be extended to practical 
models where conflict sets for a given interference model 
are known. We denote the set of SUs by S = {si, ...,sn} 
with 1 5 1 = N, and the set of channels by C = {ci, cm} 
with |C| = M. 

B. Cooperative Sensing Model 

We assume that a binary decision is made at an SU for 
each channel it senses. Let Pj(k) represent the probability 
of false alarm, i.e., the probability that a SU S{ senses 
channel k to be occupied when in fact it is idle. Similarly, 
P^(fc) represents the probability of mis-detection, i.e., 
the probability that senses channel k to be idle when 
it is actually occupied. Note that SUs outside the sensing 
range, if selected for sensing, report random sensing results. 
For instance, P^(k) = \ and P)(k) = \ if SU i is 
outside the sensing range of PU k. We assume that these 
probabilities can be learned using historical data [3], [6], 
[7]. For instance, given the location information of SUs 
and hardware parameters such as energy detection threshold 
and time bandwidth product, etc., P r l n (k) and Pj(k) can be 
calculated accordingly (see Section fV-AI for an example). 

Multi-Channel Cooperative Sensing: SUs may sense the 
licensed channels cooperatively to reduce sensing errors. To 
encourage cooperative sensing, we assume that J^. > M, 

1 Our model can be generalized to the scenario where multiple PUs access 
the same channel. 



which is common in cooperative sensing models [16], thus 
the expected number of SUs that sense a certain channel is 
at least 1, The sensing results of individual SUs are assumed 
to be independent. As mentioned earlier, due to practical 
constraints, SUs can sense a limited number of channels. 
We denote li as the maximum number of channels that SU 
Si can sense in a time slot, < h < M, for all i = 1, • • • , N 
and let l max = max^j ij. Note that li = means that 
the SU is in not in the sensing range of any channel, 
thus it cannot do any sensing and only guess the PU state 
randomly. In cooperative sensing under the multi-channel 
setting, multiple SUs choose to sense different channels and 
predict channel availability subject to the budget constraint, 
and different sensing set assignments lead to different system 
throughput across channels. We consider a centralized system 
model, where a central controller is responsible for (1) 
maintaining system parameters for PUs and SUs (2) in each 
time slot, deciding for each channel, a subset of SUs to sense 
the channel, and (3) making a global decision on channel 
availability based on the local binary decisions of SUs. Let 
Sk denote the set of SUs that cooperatively sense channel k. 
The set of all feasible channel sensing assignment policies 
are denoted by V, and defined as follows. 

Definition 2.1: Feasible assignment policy V: A set of 

sensing sets {Si, ■ ■ ■ , Sm} is a feasible assignment policy 

M 

if l| s . £ 5 fc } < li for all i, i.e., all SUs must be assigned 

fe=i 

to at most li channels to sense. 

Let Xi(k) denote the observation of channel k by SU 
Si G Sk- Further, Xi(k) = 1 represents that Si observes 
channel k to be active, while Xi(k) = represents that 
Si observes channel k to be idle. We let x(Sk) denote the 
vector of observations for channel k. Let Q = {0, 1}, and 
let Ja : — > £1 denote a general decision rule that maps 
the local observations made by a set of SUs, A C S, to 
global decision on channel activity. As the domain of /a 
will be clear from the context, we drop the subscript and 
use / instead. This decision rule applies per channel. Let 
B{k) denote the activity of channel k such that B(k) = 1 if 
channel k is occupied, and B{k) = otherwise. According 
to the definitions of false alarm and mis-detection, we define 
the conditional probability of sensing channel k to be idle 
when it is indeed idle as follows, where vector y denotes a 
particular instance of an observation vector: 

P(f(x(S k )) = 0\B(k) = 0) 

J2 P(x(S k )=y\B(k)=0), (1) 
v-f(v)=o 

where 

P(x(S k ) = y\B(k) = 0) 

n pm n (!-*/(*)). 



Similarly, we define the conditional probability of sensing 
channel k to be occupied when it is indeed occupied: 

P(f(x(S k )) = l\B(k) = 1) 
£ P(x(S k )=y\B(k) = l), (2) 



y-J(y)=i 



where 



P(x(S k )=y\B(k) = l) 

n p-^w) n p ™( fc )- 

We assume that in each time slot, a control slot T c is 
assigned for cooperative sensing, during which time a central 
controller collects P^k) and Pj(fc) from SUs, determines 
the channel sensing assignment, collects sensing results from 
SUs, and notifies an SU per channel to transmit if that 
channel is cooperatively sensed to be "idle." Note that each 
SU i only needs to send updates to the central controller 
of P^(fc), Pj(k) when their values change, e.g, when 
the location of the SU changes. Furthermore, the central 
controller only needs to compute a new assignment only 
when P^(fc), Pf(k) change. We assume T c to be a constant 
in the paper. We further assume that SUs can transmit at 
the same bit rate over each channel, and normalize this 
rate to 1. SUs are assumed to be always backlogged and 
only one of them is scheduled over channel k if sensed 
available in each time slot. Let 7To(fc) denote the probability 
that channel k is idle, which is assumed to be acquired 
accurately over time. The capacity of channel k is denoted 
by j(k) (after normalization), k = 1, ■ • • , M. We define 
= (1 - T c )7r (fc) and 2 (k) = 7 (fc)(l - Mk))- 
Following the logic in [8] and extending to the multi-channel 
case, we define the expected SU throughput over channel k 
sensed by S k . 

Ul(S k ) := (l-T c )P(B(k)=OJ(x(S k ))=0) 

= 9 1 (k)P(f(x(S k ))=0\B(k)=0) (3) 
if S k ± 0; 
Ul(S k ) := OifS fe = 0. 

where we assume that if S k = 0, no sensing is conducted 
for channel k and the channel is never accessed. Likewise, 
the expected throughput of channel k can be represented by 

Ul{S k ) := 6 2 (k)P(f(x(S k )) = l\B(k) = l) (4) 

if S k + 0; 
Ul{S k ) := e 2 (k)ifS k =%. 

Definition 2.2: System throughput: For a channel assign- 
ment {Si, • • ■ , Sm}, we define the throughput over channel 
k to be the sum of SU and PU throughput over channel 
k, denoted as U k {S k ) = Ul(S k ) + U%(S k ). The system 

M 

throughput is defined as U k (S k ). 

k=l 

Note that for a given channel sensing assignment, the 
achievable system throughput is determined by the decision 



rule /. In this paper, we apply the optimal Bayesian deci- 
sion rule proposed in [8] to each channel respectively, to 
obtain the optimal expected system throughput. Formally, 
for each channel k and an observation vector y by S k , 
if e 2 (k)P(x(S k ) = y\B(k) = 1) > 0i(AOP(sc(S fc ) = 
y\B(k) = 0), the decision on channel k is "occupied", and 
the contribution to throughput is 62(k)P(x(S k ) = y\B(k) = 
1) ; otherwise, the decision on channel k is "idle" and the 
contribution is 9i(k)P(x(S k ) = y\B(k) = 0). 

C. Problem Formulation 

We formulate the optimization problem to maximize the 
system throughput, including PUs and SUs on all channels, 
as follows: 

M 

Problem (A): max V U k (S k ), 
{S 1 ,-,S M }ePfe=i 
where the Bayesian decision rule is implicit in the definition 

of U k (-). 

Our goal is to decide the optimal channel sensing assign- 
ment to maximize system throughput. We adopt a common 
assumption that PUs can tolerate interference to a certain 
extent, which may appear in the form of a constraint as in [3], 
[9] and our earlier paper [8] for the single channel setting. 
In the future, we plan to extend our solution presented in 
this paper to Problem (A) with explicit constraints on PU 
throughput. 

We assume that the system is static and the optimization is 
done in a single time slot. Note that the solution of the static 
assignment would apply to multiple time slots if P? n (k) and 
Pf{k) do not change over time, or if changes occur over a 
much slower time scale. 

III. Hardness of the Problem 

In this section, we will show that Problem (A) is strongly 
NP-hard [14], by a reduction from Product Partition, which 
is NP-complete in the strong sense [1]. The Production 
Partition problem is defined as follows: Given N positive 
integers oi, a 2 , on, is there a subset X C J\f := 
{1,2, N} such that f] a » = II a 4 ? 

We reduce Product Partition to the following subproblem 
of Problem (A), with M = 2, Pj(l) = P|(2) = for all 
i, P^(l) = P4(2) := Pr n for all i, and k = 1 for all i, 
7 (1) = 7 (2) := 7 , 7r (l) - tt (2) := tt , (1 - T c )ir := U 
7(1 - 7r ) := 9 2 , and 61 = 2 . 

Let (Si, S 2 ) denote a solution to this subproblem. Without 
loss of optimality, we can assume Si and S 2 form a partition 
of the set of SUs, i.e., Si US 2 = S and Si HS 2 = 0. The ex- 
pected system throughput can then be easily determined us- 
ing the Bayesian rule as Ui(Si) = 6i+6 2 {l- f[ P^J and 



U 2 (S 2 ) = 9 1 - 
2f?i +( 



2(1- 



ieSi 



[ P^j). Problem (A) then becomes: 

2(2- ( n pl+ n p, 

further equivalent to min ( ] P^ + 



max 

SiCS 



which is 
PL) since 



2#i + 2#2 is a constant. We then establish the strong NP- 
hardness of Problem (A) by showing that this new problem 
is strongly NP-hard. 

Proposition 3.1: Problem (A) is strongly NP-hard. 

Proof: By the above argument, it suffices to prove that 
the subproblem, min ( ] [ P 7 l n + P, l n )> is strongly 

Sl ^ s s^Si Si es\Si 
NP-hard. Given an instance of Product Partition with pa- 
rameters oi," - ,ajy, we reduce it to an instance of this 
subproblem as follows: let P r l n = ai/HY, i = 1, • • • ,N, 
where r is the smallest integer such that P l m < 1 for all i = 
1,...,N. This reduction can clearly be done in polynomial 
time. Furthermore, if the re is a subset X C J\f, such that 

\ di = Oi = / Y[ a u th en the optimal solution 

iex iej\f\x y ieM 

to the subproblem is 2 / Y[ Pm> an ^ vice-versa. Hence if 

there is polynomial time algorithm to the subproblem, the 
Product Partition problem can be determined in polynomial 
time as well, which contradicts the fact that Product Partition 
is strongly NP-complete. ■ 
Since Problem (A) is strongly NP-hard, no pseudo- 
polynomial time algorithms exist unless P = NP [14]. We 
will propose a matching-based approximation algorithm that 
has theoretical lower bound in Section [TV] 

IV. Approximate Solutions 

In this section, we propose an efficient approximation 
algorithm for Problem (A). We first prove an upper and 
a lower bound on system throughput. We then propose 
a matching-based approximation algorithm. By exploiting 
the structural properties of the problem and the bounds on 
system throughput, we show that the algorithm achieves an 
approximation ratio of at least i/i, where /i <E [1,2] is a 
system parameter and will be defined later. 

A. Property of the System Throughput 

We will show the range of the system throughput Uk(-) 
in the following lemma. 

Lemma 4.1: For any SU s 2 and channel we have 

e 1 (k) + e 2 (k) > U k (m) > max{e 1 (k),6 2 (k)}. 
Proof: 

Uk({ Si }) = r a ^[e 1 {k){l~P}{k)),e 2 {k)P l m {k)} 

+ max { 0i {k)P} (k) , 9 2 (k) (1 - P^ {k)) } 

> max{<9i(fc),<9 2 (fc) j 

Furthermore, it is clear from the definition of Uk{-) that 
6»i(fc) + 2 {k) > U k (si), since at most both PU and SU can 
achieve their full capacity. ■ 

B. A Matching-Based Approximation Algorithm 

In this section, we propose a maximum weighted matching 
(MWM) [15] based algorithm to Problem (A). We first 
provided a detailed description of our algorithm (see Algo- 
rithm Q]), and then establish its approximation factor. 

The algorithm starts with constructing a complete and 
weighted bipartite graph (lines 2-4), where for each channel 



Algorithm 1 A maximum weighted matching based algo- 
rithm for maximizing the system throughput across channels 
Input: N, M, T c , Tv {k), j(k) for all k; h for all i; P^(fc), P}(k) 
for all i and k 
Output: U and S k for all k 
1: S k <- for all k 

2: V<-{s{,--- ,s[ 1 ,--- ,s]v,-- - ,s^}U{ci,-- - ,c u } 
3: /' < U. :.. ..v *=!,...,* {UjU (»*.<*)} 
4: G <- (V, E) 

5: w(4,c k ) <h- U k ({Si}), Vi = 1, • • ■ ,N,j — 1, • • ■ ,U, k = 
1,.-. ,M 

6: M <— a maximum weight matching in G 
1: S k <- { Si : (si,c k )€M},Vk 
8: R {si : si is not matched in M} 
9: for all sj € R do 

10: fc'^arg max \u k (S k U {st}) - Uk(Sk)] 

11: Sk* <~ S k * U {si} 

M 

12: U<- E Uk(Sk) 

k = l 

13: Ui •(- max£4(S) 

14: if Ui > U then 

15: f7 Ui 

16: fc* argmax Uk(S) 

17: S, Sk\- Vfe ^ k* 



k, a vertex is constructed, and for each SU Sj, ij vertices 
are constructed corresponding to the li copies of the SU, 
denoted as s\, j = and for any pair of vertices s\ 

and Cfc, there is an edge connecting them. The weight of an 
edge (s 3 i7 Ck) is then defined as w(s\,Ck) = Uk({si}) (line 
5). 

A maximum weight matching in the bipartite graph is 
then found (line 6), and for each edge (s^,Cfc) in the 
matching, SU Sj is assigned to sense channel Cfc. A greedy 
heuristic is applied for determining the assignment of the 
remaining copies of SUs to channels (lines 8-12). Basically, 
the remaining copies are first sorted in an arbitrary order, 
and a copy of is assigned to the channel that provides the 
maximum marginal improvement of the system throughput 
among all the channels not assigned to s; yet. This scheme 
is then compared with another scheme for which all SUs are 
assigned to a single channel that gives maximum throughput 
(line 13). The algorithm outputs whichever scheme provides 
a larger system throughput. 

We then analyze the complexity of Algorithm Q] which is 
dominated by computing the maximum weighted matching 
and evaluating the throughput function Uk(-). It is shown 
in [8] that for a given sensing set Sk, Uk(Sk) can be eval- 
uated using a dynamic programming algorithm in pseudo- 
polynomial time. Let Q denote the time complexity for one 
evaluation of Uk{-)- Note that the total number of such eval- 
uations is bounded by Nl max M. Therefore, the time com- 
plexity of Algorithm |T] is 0(Nl max MQ + (Nl max + M) 3 ). 

To establish the approximation ratio of Algorithm Q] 
we first construct a maximal matching called M_Gdy that 
approximates the MWM and captures two key aspects: 1) 



SUs may have different sensing abilities for each channel; 2) 
channels are competing for SUs with limited sensing budget. 
We first prove a lower bound on the system throughput using 
the sensing assignments determined by M_Gdy, which is then 
used to prove the approximation ratio of Algorithm [TJ The 
matching is constructed as follows: 1) Partition the channel 
set C into groups indexed by SU, and each group is labeled 
as Ci that includes all channels k with Uk({si}) > Uk({sj}) 
where j ^ i. Ties are randomly broken. Let denote the 
size of Cj. 2) Sort the channels k in each group Cj by U% 
in descending order, where U% = min ie s J7fc({sj}). 3) Pick 
the first li channels from each group Cj (the set is labeled as 
Cj 4 ) and assign SU i to sense these channels. 4) Randomly 
assign an unused SU copy to each of the rest channels. We 
will next show a lower bound on the system throughput 
using M_Gdy in Lemma |4~2l We define A; = min{Zj, r^jri, 
Pi = min fe c i t -jjfc where E/£ = maxj Uk({si}). Note that 
by Lemma l44l pi G [1,2] for all i. We have the following 
performance bound, where p = 1 + mirijgs Xi{pi — 1), and 
\M_Gdy\ denotes the system throughput using the sensing 
assignments determined by M_Gdy. 

Lemma 4.2: M_Gdy achieves a system throughput no less 
than/x£ fe ^°- 



Proof: 
\M.Gdy\ 



Ei £ s [E 



Ei G S [E 
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EiGS Efc 



ec, 
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E l6S [(p;-i)E fceC ^ ] 



> i + 
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EiGS EfcGC 



ec, 



Eigs E 



GS ^*6C, 

1 + min Ai(pi — 1) = p 

iGS 



Based on Lemmas l4.1l and l4T2l we show the approximation 
ratio of Algorithm [Tj in Proposition 14.31 

Proposition 4.3: Algorithm [TJ achieves at least a fraction 
of 7}p of the optimal system throughput for Problem (A). 

Proof: Let OPT be the optimal solution, and ALG be 
the solution by Algorithm [TJ to Problem (A). By Lemma |4~T1 
we know that 



Eh uj > E/t 

OPT ~ 



max{6>i(fc),6> 2 (fc)} 1 
£k *i (*)+&(*) " 2' 



(5) 



Since ylLG is an outcome at least as good as maximum 
weight matching and M_Gdy is a matching we construct in 
a greedy way, we have ALG > \MjGdy\. By Lemma l4~2l 
we can achieve ^\l i - ■ 

Remark 1: We note that when 9i{k) > 9 2 {k) or 9 2 (k) > 
6\{k), we can achieve a solution close to the optimal by 
Algorithm [TJ since Equation (0 becomes close to 1. Only 
when 6i(k) and 62(h) for all k are close, Equation (0 is 
only right above |. Also, if SU's sensing abilities across 



channels vary in a large range, or the sensing budgets of 
SUs are large, the ratio will be close to 1 since pi, A, will 
be large, respectively. 

Remark 2: In the proof of Proposition l4.3l we have ignored 
the greedy heuristic applied to the copies of SUs not included 
in the matching. Hence the result established above only 
provides a lower bound on the performance of our algorithm. 
Proving a tighter bound for the algorithm that incorporates 
the greedy heuristic is part of our future work. 

V. Simulations 

In this section, we study the performance of our algorithm 
through simulations by comparing Algorithm [TJ (MWM) 
with a random sensing assignment algorithm, and a greedy 
algorithm (defined next). In the random algorithm, the copies 
of SUs are randomly assigned to PUs. The greedy algorithm 
works as follows: for each PU k, the set of SUs are first 
sorted by -P^(fc) + Pj(k) in a non-decreasing order as its 
preference list. In each round, a random permutation of the 
set of PUs is applied. The algorithm then goes through the 
PU list, and for each PU k, a copy of the SU, say s,;, with 
the lowest P? n (k) +Pj(k) among the remaining SUs, which 
has not been assigned to k before and has remaining copies, 
is assigned to k. Repeat this procedure till all copies of SUs 
have been assigned. 

A. Simulation Setting 

The following parameters are fixed throughout the simu- 
lations. We consider a 100 x 100 area, where the locations of 
M PUs are randomly generated. For each PU k, its maximum 
power level is randomly chosen between 1 and 10, and 7To(fc) 
are randomly generated in [0, 1]. We also set T c = 0.2 fixed. 
In each of the 100 runs of the simulation, we apply the model 
proposed in [11] to generate i^(fc) and Pj(k). The details 
are in our online technical report [17]. 

B. Simulation Results 

The simulation results are shown in Figure [TJ Note that 
we do not restrict Yli h ^ M in our simulations. If PU k 
is not assigned any SU for sensing (St = 0), the system 
throughput on channel k is 62(h) (Definition 12.2b . In all the 



di(k) + e 2 (k) 



as the upper bound for 



figures, we plot J2k 
the optimal solution. 

In Figure [T(a)l we fix M = 20, l max = 3, and vary N from 
4 to 20. For each PU k, j(k) in generated randomly in [1, 3] 
and then fixed over all 100 runs. We choose this range since 
the average PU throughput is usually larger than the unit SU 
throughput. For each SU i, k is randomly generated between 
1 and lmax and fixed over all the runs. The simulations results 
are averaged over all 100 runs. We observe that AlgorithmQ] 
achieves significant improvement over the random and the 
greedy algorithms for all N, although the gap shrinks as N 
increases. For instance, the system throughput of AlgorithmQ] 
is 24% larger than that of the greedy algorithm when = 4 
and it decreases to 16% when N = 20. When more SUs 
join the network, the random and the greedy algorithms have 
more chance to choose "good" SUs. The greedy algorithm 
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Fig. 1. System throughput achieved by our algorithm, greedy algorithm and random algorithm 
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is comparable to the random algorithm when N is small. 
However, it wins over the latter when N > 12. This indicates 
that the sorting step in the greedy algorithm helps PUs pick 
the "right" SUs, which is more useful when N is large. 
Besides, the performance of Algorithm Q] reaches 95% of 
the upper bound of the optimal solution when N = 20. 

In Figure [1(b)] AI, N are fixed to be 20 and 8, respectively, 
and we vary l max from 1 to 5. 7(fc) is again generated 
randomly in [1,3] and fixed over all 100 runs. Similar to 
Figure |l(a)| Algorithm [T] outperforms both the random and 
the greedy algorithms, and the greedy algorithm outperforms 
the random algorithm when l max > 3. When l max = 4, the 
system throughput of Algorithm Q] is 25% better than that of 
the greedy algorithm, which is the largest gap in the figure. 
An interesting observation is that the expected number of 
SU copies when l max = 4 (N = 8) is equal to M = 20, 
thus every PU is assigned an SU on average. When there 
is more supply (SUs) than demand (PUs) or more demand 
than supply, the performance gap between Algorithm Q] and 
greedy algorithm, random algorithm is not so significant. 

In Figure [1(c)] we fix M = 20, N = 8, l max = 3, and vary 
the range of the channel capacity j(k). For instance, [1,2] 
means all channel capacities are randomly generated between 
1 and 2. Algorithm Q] is constantly better than the other two 
algorithms. The gap first increases as the channel capacity 
increases (from 18% to 34%) till -f(k) S [1, 2], and decreases 
thereafter (7% at j(k) e [1,5]). When the channel capacity 
is comparable to unit SU capacity, the choice of SUs for 
sensing does not affect the system throughput significantly; 
When the channel capacity dominates the system throughput, 
the choice of SUs again loses its leading role. Thus the 
largest gap appears in the middle. 

VI. Conclusion 

In this paper, we investigate the problem of throughput 
maximization using cooperative sensing in multi-channel 
CRNs, where each SU can only sense a limited number of 
channels with various sensing capabilities, due to time or 
energy constraints. We show that under the optimal Bayesian 
decision rule, the channel sensing assignment problem is 
strongly NP-hard. A matching based algorithm is then pro- 
posed with an approximation ratio that is at least ifx where 
fj, G [1,2] is a system parameter. Our numerical results 
demonstrate that our algorithm performs significantly better 
than the a random channel sensing assignment algorithm and 
a greedy algorithm. As part of our future work, we plan 



to establish a tighter performance bound for our algorithm 
enhanced with a greedy heuristic, and consider the system 
throughput maximization problem with extra constraints on 
the PU throughput. 
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